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Abstract 

Einstein's equations in a tetrad formulation are derived from a lin- 
ear theory in flat spacetime with an asymmetric potential using free field 
gauge invariance, local Lorentz invariance and universal coupling. The 
gravitational potential can be either covariant or contravariant and of 
almost any density weight. These results are adapted to produce uni- 
versally coupled massive variants of Einstein's equations, yielding two 
one-parameter families of distinct theories with spin 2 and spin 0. The 
theories derived, upon fixing the local Lorentz gauge freedom, are seen 
to be a subset of those found by Ogievetsky and Polubarinov some time 
ago using a spin limitation principle. In view of the stability question for 
massive gravities, the proven non-necessity of positive energy for stabil- 
ity in applied mathematics in some contexts is recalled. Massive tetrad 
gravities permit the mass of the spin to be heavier than that of the spin 
2, as well as lighter than or equal to it, and so provide phenomenological 
flexibility that might be of astrophysical or cosmological use. 
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1 Introduction 



The project of deriving a relativistic gravitational theory using considerations 
such as an analogy to Maxwellian electromagnetism, the universal coupling of 
the gravitational field to a combined gravity-matter energy-momentum complex, 
and the requirement that the gravitational field equations alone (without the 
matter equations) entail energy-momentum conservation was a part of Einstein's 
search for an adequate theory of gravity in 1913-15 [1-4]. This history was 
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subsequently downplayed by Einstein [4,5], and the above ideas later came 
to be associated with the (allegedly non-Einsteinian) field-theoretic approach 
to gravitation more commonly associated with particle physics, where it was 
brought to successful completion in the 1950s [6,7]. 

For some time it was believed that there was a unique massive generalization 
of General Relativity satisfying universal coupling [8,9]. This theory has the 
same mass for the spin as for the spin 2, as well as a specific nonlinear algebraic 
self-interaction such that the DeDonder-Fock harmonic condition is satisfied by 
the effective curved metric when Cartesian coordinates are used for the flat 
background metric. 

A few years ago it was found that there are, at least, two one-parameter fam- 
ilies of universally coupled massive theories [10] . These involve either a covariant 
or contravariant symmetric tensor potential of almost any density weight; the 
resulting theories involve mass terms conveniently presented in terms of cor- 
respondingly densitized curved metric or inverse metric. All of these theories 
are included in the 2-parameter family of theories found by Ogievetsky and 
Polubarinov in terms of their spin limitation principle [11,12], which imposed 
axiomatically an auxiliary condition to eliminate the spin 1 and one spin 
degrees of freedom. 

An important but often neglected aspect of deriving Einstein's equations 
as the field equations for a self-interacting massless spin 2 field theory is the 
inclusion of spin i matter, that is, fermions. This problem was addressed using 
in effect an orthonormal tetrad formalism in an unjustly neglected paper by 
Shirafuji [13]. Going back to Weyl in 1929 [14,15], it is often claimed that 
coupling of spinors to the curved space-times implied by General Relativity 
requires the introduction of an orthonormal tetrad. While that claim was proven 
false in the mid-1960s [11, 16, 17], nonetheless in some contexts it is convenient 
to represent the gravitational potential using an orthonormal tetrad, especially 
(but not only) if spinors are present, in order to avoid technically demanding 
nonlinearities. The interesting question of including spinor matter is set aside for 
another occasion, however, in favor of deriving mass terms that can accompany 
the Einstein tensor. 

Using a tetrad formalism, I will exhibit two one-parameter families of mas- 
sive gravities which have not previously been derived from universal coupling. 
Einstein's equations have been derived using universal coupling with a tetrad 
and an independent connection [18, 19]. While this paper draws on themes from 
those papers, it aims at deriving Einstein's theory and then massive generaliza- 
tions, treating the connection implicitly as some function of the tetrad. It seems 
to me that the present derivation, which does not aim to arrive at supergravity 
or even gravity with spinors, has the advantage of requiring little knowledge 
of the answer beforehand, as well as the simplicity of fewer field variables. It 
is therefore quite useful in the context of discovery, not just the context of 
justification. 
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2 Stability or Empirical Adequacy Questions for 
Massive Gravity 

Prom the early 1970s until roughly 2000 (as dark energy was becoming ac- 
cepted), the conventional wisdom about (Lorcntz-invariant) massive variants of 
General Relativity was that they came in two varieties, which succumb to dif- 
ferent problems [9,20-23]. The Pauli-Fierz pure spin 2 theory has 5 degrees of 
freedom per spatial point (at least to lowest order), avoiding a wrong-sign spin 
field by making it infinitely massive. But the van Dam-Veltman-Zakharov 
discontinuity shows (at least perturbatively) an observable discontinuity in the 
massless limit, by which pure spin 2 massive gravities already were experimen- 
tally excluded: light spin 2 graviton theories do not approximate Einstein's 
equations. The other possibility, involving spin 2 and a wrong-sign spin at 
every point, is not expected to have such a discontinuity, but instead is believed 
to be violently unstable. The desire to avoid negative-energy degrees of freedom 
motivates gauge invariance for linear theories [18, 24] and thus leads to massless 
spin 2, that is, Einstein's equations. 

In the last decade or so, however, the consensus has weakened in both cases, 
partly because dark energy has provided empirical evidence that Einstein's equa- 
tions might be wrong on long distance scales and partly because explanations 
of dark energy tend to involve weakening assumptions of energy positivity. One 
finds claims and counterclaims regarding the non-perturbative smoothness of 
the massless limit of Pauli-Fierz pure spin 2 theories [25-34] . The infinite-mass 
cases below are of this sort. One also finds both increasingly frequent favor- 
able mentions of spin 2-spin theories [35-40] and, sometimes independently, 
increasing doubts that negative energy degrees of freedom are always bad. Two 
considerations warrant giving spin 2-spin theories enough consideration to 
assess by detailed calculation. First, the old consensus against such theories 
appears to be motivated partly by a false converse of a true theorem. While 
positive energy is sufficient for stability in mechanics, demonstrably it is not nec- 
essary, as is known in nonlinear stability theory and plasma and fluid physics 
[41-45]. Much depends on the resonances and nonlinear interaction terms, at 
least for classical finite-dimensional theories and field theories in a box (which 
seem to be the state of the mathematical art), and presumably on the disper- 
sion relation(s) as well. Of course things might well be harder in field theory 
(in infinite space), and harder still under quantization. It is noteworthy that 
one also sees Kolmogorov-Arnol'd-Moser-type theorems being proved for partial 
differential equations [46,47]. But clearly negative energy degrees of freedom 
do at least make the question of stability significant. The second consideration 
is that there are (previously unnoticed, evidently) reasons to wonder whether 
there really are negative energy degrees of freedom at high wavenumbers. While 
the claim that the spin has negative energy goes back to Pauli and Fierz, the 
calculation turns out to be nontrivial. For the Hamiltonian in terms of the true 
degrees of freedom, taking the simplest cases with quadratic dependence on the 
shift vector [9, 10,48], eliminating the shift vector using its equation of motion 
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gives a positive cross-term ^ TL^ mixing the wrong-sign spin and right-sign 
heUcity part of the spin 2. The reduced Hamiltonian density contains not 
only derivatives of the canonical momenta due to eliminating the shift, but also 
higher spatial derivatives due to eliminating the lapse, and so docs not fit the 
usual form. (The Lagrangian density, by contrast, becomes formally nonlocal, 
much as in ([48]), and hence not clearly simpler.) Apart from detailed calcula- 
tion, it isn't just obvious that there really arc negative energy degrees of freedom 
at high wavenumbers, where the strange new terms will be important. Negative 
energies for low wavelengths only, it has been suggested for other theories, might 
permit mctastability [49-51]. There are reasons to worry about negative energy 
problems, but not such compelling reasons as to obviate detailed investigation. 
The theories derived here include both pure spin 2 and spin 2-spin examples, 
so resolving the objections to either sort of theory would vindicate some of these 
theories. 

3 Curvature from Universal Coupling and Gauge 
Invariance: Covector Density Potential 

One can consider what theories result from a covector potential (with an ad- 
ditional local Lorentz index) with free field gauge freedom and coupling to the 
total stress energy tensor. Any density weight will be permitted, except for a 
singular case. The massless case will be treated first, yielding Einstein's equa- 
tions. Later a mass term will be introduced. 

3.1 Free Field Action 

For the massless theories, one assumes an initial infinitesimal invariance (up 
to a boundary term) of the free gravitational action, along with local Lorentz 
freedom. For the later derivation of massive theories, the gauge freedom will be 
broken by a term algebraic in the fields, but the derivative terms will retain the 
gauge invariance. The entire theory will retain local Lorentz freedom. 

Let Sf he the action for a free covector tensor density 7^ (of density weight 
—w, w j) in a. spacetime with a flat metric tensor 77^,^ in arbitrary coordinates 
in four space-time dimensions. The flat metric tensor r/y^i/ is equivalent as a 
geometric object to the weight — 2w tensor density ^^j^; the two are related by 

flf.. = ri,.v^-^'"- (1) 

One notes that the forbidden case w y^ j makes 77^;/ non-invertible as a function 
of the metric: 77^,/ -^/^ ^^^^ determines only the null cone, not a full metric 
tensor. The weight — 2w metric density can be built from the weight —w 
cotetrad density fj^ by the relation 

Vtiu = V^VabV^, (2) 
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where rjAB is just the matrix diag{—l, 1, 1, 1). fj^ will represent the geometry 



of space-time in this derivation. One can therefore write the action for the free 



gravitational field as <S'/[7f , r)f^]. This action also involves the constant matrix 



VAB- 

The torsion- free metric-compatible covariant derivative is denoted by d^, so 
daVnv = and dafliiv = 0. One recalls that covariant and Lie differentiation of 

densities involves an extra term due to density weight. A (1,1) density of 
weight V is representative. The Lie derivative is given by [52] 

£i4>"p = e^h^^ -^^^"'M +^X'0 +^<^^c^M • (3) 

The ry-covariant derivative is given by [52] 

d^4>1 = 4>U +^?r«^ - Cr?^ - v4>'iTi^. (4) 

Here F^^ are the Christoffel symbols for r]nv Once the curved metric g^v is 

defined below, the analogous g-covariant derivative V with Christoffel symbols 
{^^} follows. On account of the local Lorentz freedom, it is not the case that 
daf]'^ vanishes. While one could define another covariant derivative that corrects 
this problem, such a construction will not be needed here. For the dcnsitized 
flat metric 77^1/, the oppositely densitized inverse flat metric is f}'^" . Likewise, the 
inverse of the densitized curved metric 5^1/ to be defined below is g^" . 

The potential problem of wrong-sign degrees of freedom can be avoided for 
massless theories using suitable gauge freedom. Let us require that the free 
field action 5/ change only by a boundary term under the infinitesimal gauge 
transformation 7^ — ^ 7^ + ^7^, where 

One anticipates that a connection between w and c will emerge. This expres- 
sion does not include the possible term S," dufj'!} , which would spoil local Lorentz 
invariance of this gauge transformation and also the form of the resulting gener- 
alized Bianchi identity as the divergence of an expression linear in |||- . Neither 

does it include Vi^pfigri^^di,^^^ , which does not appear to be helpful for arriving 
eventually at an expression of the form (5.g^,y ^ 

For any Sj (quasi-)invariant in this sense under (5), a certain linear combi- 
nation the free field equations is identically divergenceless. The action changes 
by 



5Sf = J dt^x 



(6) 



plus an irrelevant boundary term. Integrating by parts, letting have compact 
support to annihilate the boundary terms (as we shall do throughout the paper), 
and making use of the arbitrariness of one obtains the identity 



-8.(.5^^ + <^Hi|^l=o. (T) 
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This is the generalized Bianchi identity for the free theory on account of the 
gauge freedom. 

Besides this identity, there is also the local Lorentz invariance of the free 
theory (not to mention the full theory) to consider. The free field action Sf, 
under an infinitesimal local Lorentz transformation 

= ^iv^, (8) 



is left unchanged. The matrix field fi^ is antisymmetric when an index is moved 
by rjAB or its inverse. The resulting generalized Bianchi identity is 

SSf _ SSf ^ 
l^.[B—T]+V^.lB—i^=0. (9) 



3.2 Cotetrad Stress-Energy Tensor Density 

If the energy-momentum tensor is to be the source for the gravitational poten- 
tial 7^, consistency requires that the total energy- momentum tensor be used, 
including gravitational energy- momentum, not merely non-gravitational ( "mat- 
ter") energy- momentum, for only the total energy-momentum tensor is diver- 
genceless in the sense of di, [53,54], or, equivalently, in the sense of a Carte- 
sian coordinate divergence. To obtain a global conservation law, one needs a 
vanishing coordinate divergence for the 4-current. Thus a vanishing covariant 
divergence (in terms of a curved connection) for a two-index energy-momentum 
complex is not useful. 

An expression for the total energy-momentum tensor density can be derived 
from S using a tetrad (or cotetrad, or densitized tetrad, or densitized cotetrad) 
that corresponds to a fiat metric tensor 77^,^; for brevity, the entity will be called 
a flat (co)tetrad, even if it is densitized. One has = fj^riABflu , where 
Vab = diag{—l, 1, 1, 1). Let the action S depend on the flat cotetrad fj^, the 
(asymmetric) gravitational potential 7^, and matter flelds u. Here u represents 
an arbitrary collection of dynamical geometric objects (in the sense of having a 
set of components at each point in each local chart and a transformation rule 
in the overlap regions [55,56]), having perhaps some coordinate indices but no 
Lorentz indices. Under an arbitrary inflnitesimal change of coordinates with 
compact support (so that boundary terms can be discarded without notice), 
described by a vector fleld the action changes by the amount 

with boundary term BT vanishing due to compact support of But S* is a 
scalar, so SS = 0. Integrating by parts, discarding more vanishing boundary 

terms, and using the arbitrariness of the vector field gives a generalized 
Bianchi identity, from which the conservation of stress-energy can be inferred. 
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To have a moderately nontrivial sort of matter for dcfiniteness, let u be a single 
contravariant vector field u^. The resulting identity is 



(11) 

The first term jM^dafl^ is neither part of a divergence nor vanishing on-shell 

(that is, using the matter and gravitational equations of motion); it also violates 

local Lorentz invariance. Thus it will need further attention. 

For the local Lorentz freedom, a similar derivation to that for Sf above 
shows that 

SS ^ 5S ^ 

which implies that the antisymmetric part of the stress-energy tensor vanishes 
when gravity is on-shell. Using the local Lorentz identity in the coordinate 
identity and rearranging some terms gives 



ss ^ _A { ~A ss Tn~A SS \ ss ^ ,, „ 



1 = 0. 



Letting matter u and gravity 7^ satisfy their equations of motion gives a con- 
served symmetric stress-energy tensor 

This quantity is an energy-momentum tensor density for matter and gravita- 
tional fields. As usual with Rosenfeld-style derivations using a derivative with 
respect to a flat metric, one relaxes flatness while taking the functional deriva- 
tive and then restores flatness later [57,58]. 



3.3 Full Universally-Coupled Action 

We seek an action S obeying the plausible physical postulate that (invertible 
linear combinations of) the Euler-Lagrange equations be just (invertible linear 
combinations of) the free field equations for 5/ augmented by the total energy- 
momentum tensor. A simple way to impose this requirement using a densitized 
covector to represent the flat space-time metric and another one to represent 
the gravitational potential is: 

SS__5S^_X5S_ 
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where A = — v327rG. The expression on the right involves a trace- altered rel- 
ative of the stress-energy tensor. (One could shift the trace and obtain the 
stress-energy tensor on the right side and a trace-altered relative of the Euler- 
Lagrange equations if one wished.) One is free to make a change of variables 
in S from 7^ and f]^ to the 'bimetric' (really "bitetrad," but that word is not 
catchy) variables and fj-^, where 

9t = vt- ^7^ (16) 

(From g-^ one can then define the metric gf^y by matrix algebra and then define 
the g-covariant derivative V in the usual way. We shall have little need for 
explicit use of V, however.) Equating coefficients of the variations using the old 
variables and the new variables in terms of the old gives 



5S _ SS ^ 6S 



aI'T = J^\~9 + (17) 



for Sf/-^ and 



SS _ A SS 



(18) 



for ^7^. Putting these two results together gives 

SS ^SS 



^7-^l7 = A^|ff-2-^ (19) 



which splits the stress-energy tensor into one piece that vanishes when gravity is 
on-shell and one piece that does not. Using this result in the universal coupling 
postulate gives 

A SS SSf 

Multiplying by the flat tetrad, taking the divergence, and recalling the free 
fleld theory's Bianchi identity from gauge invariance gives 



9AtT^\~9 + cS>:;f,t^\~g\=Q. (21) 



Above c was a free parameter, so it now makes sense, in view of the form of 
the stress-energy tensor, to set c = —w. By playing with the density weights 
to account for the trace term, one can show that the resulting equation is, in 
terms of the non- weighted tetrad 77^ (with no"), 

0.{v^^\9)=0. (22) 
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Raising an index gives 



9Jv''^t-t\~9]=0. (23) 



SrjA 

There being no spinors present in the theory, one can choose for the ciirvcd 
tetrad to appear only through the curved metric g^i,, as one knows occurs 
in Einstein's equations. Local Lorentz invariance does not exclude dependence 
on, for example, g^fj'X- Such dependence will appear below in the mass term 
as well. But it is difRcult to envision how a term that depends on both the 
curved and flat tetrads could satisfy the previous equation as an identity. So 
restricting the action to depend on the curved tetrad only through the curved 
metric might not be an additional assumption after all. With dependence on 
the curved tetrad restricted to dependence through the curved metric, it follows 
as a theorem that dependence on the flat tetrad is also only through dependence 
on the corresponding flat metric 77^^ [59] . One then has 

One might wonder what has been gained using the tetrad formalism instead of 
a metric formalism. (A tempting wrong answer is that one acquires the ability 
to have spinors in the theory. But spinors were already admissible at the cost 
of technical complexity and unfamiliarity, as noted above.) The benefit will 
appear when the mass term is introduced, because the tetrad or cotetrad will 
permit a simple derivation of a mass term. The density weight of the tetrad 
and the type (contravariant or covariant) of tetrad will determine the specific 
mass term obtained. 

The flat metric appears, after some changes of variables, only in a (sym- 
metric) expression with identically vanishing divergence. This fact, combined 
with the split of the stress-energy tensor into two parts, shows that the gravi- 
tational field equations alone entail conservation of energy-momentum, without 
any separate postulation of the matter equations. The quantity -^^\g, being 
symmetrical and having identically vanishing divergence on either index, neces- 
sarily has the form 



5S 



\g = ^dpd^ (^A^MKl + TV^MM^ + 3^7]^"" (25) 



[60] (pp. 89, 429), where M'^p'^" is a tensor density of weight 1 and i? is a con- 
stant. This result follows from the converse of Poincare's lemma in Minkowski 
spacetime. M^'"^'^ cannot be chosen arbitrarily, but rather must be chosen so 
that the term 4^ is accommodated. 

Gathering all dependence on 77^^ (with g^^, or equivalcntly g^^, independent) 
into one term yields S = <S'i[^^y,u] -|- S'2[g^i/, jy^^i/, w]. Using the effective curved 
metric density ^y^^, one can define an effective curved metric by 5^1/ — g^iv\/—g~^ 
(where I = 2w) and an inverse curved metric density g'^" . For 6*1, we choose 
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the Hilbcrt action for general relativity plus minimally coupled matter and a 
cosmological constant: 



As is well-known, the Hilbert action is the simplest (scalar) action that can be 
constructed using only the metric tensor. (One could also admit terms with 
higher derivatives of the curved metric.) If the gravitational field vanishes ev- 
erywhere, then the gravitational action ought to vanish also. In the massless 
case, the result is that B = A/lGnG. For the massive generalization below, 
the gauge-breaking part of the mass term will introduce another zcroth order 
contribution that also needs canceling. One could also let the matter couple to 
the Riemann tensor for g/^i, or allow higher powers of the Riemann tensor into 
the gravitational action, if one wished. For the massless case, one might set 
A = [8] . One easily verifies that if 



then j^lg has just the desired form, while 52 does not affect the Euler- 

Lagrange equations because = and = identically [6] . The coefficient 
B of the 4- volume term is naturally chosen to cancel any other zeroth order term 
(such as from a cosmological constant) in the action, so that the action vanishes 
when there is no gravitational field. The boundary term is at our disposal; if 
a'^ is a weight 1 vector density, then 5 is a coordinate scalar. It is doubtful 
that the unmodified Hilbert action is best, given its badly behaved conservation 
laws with the factor of 2 problem [61,62]. In summary, the universally coupled 
action for the massless cotetrad density is 



for which the first few terms give Einstein's equations and the last few give no 
contribution to the field equations. 

4 Massive Cotetrad Density Theories 

The present goal is to generalize the derivation above to yield one or more mas- 
sive finite-range variants of Einstein's equations. Such field equations would 






(27) 




(28) 
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relate to Einstein's in much the way that Proca's massive electromagnetic field 
equations relate to Maxwell's [63]. There are some differences, however, includ- 
ing the lack of uniqueness for massive gravity, even given the plausible postulate 
of universal coupling. It will turn out that every real number (with one excep- 
tional forbidden case) gives a distinct massive variant of General Relativity. 

One expects that the mass term for a free field be quadratic in the potential 
and lack derivatives. The free field action Sf is now assumed to have two 
parts: a (mostly kinetic) part 5/0 that is invariant under the erstwhile gauge 
transformations as in the massless case above, and an algebraic mass term Sfm 
that is quadratic and breaks the gauge symmetry. We seek a full universally 
coupled theory with an action S that has two corresponding parts. The two 
parts of 5 = So + Sms are the familiar part 5o (yielding the Einstein tensor, the 
matter action, a cosmological constant, and a zeroth order 4-volumc term) and 
the new gauge-breaking part Sms which also has another zeroth order 4- volume 
term. As it turns out, the mass term is built out of both of the algebraic part 
of Sq (the cosmological constant and 4-volume term) and the purely algebraic 
term Sms- Moreover, every part of the action is locally Lorentz-invariant. 

Requiring Sfo to change only by a boundary term under the variation S-y^ = 
Vv^n^'^ — wfj'^dv^," implies the identity 

Again we postulate universal coupling in the form 

5^_5S^_\5S_ 

Changing to the bimetric variables implies, as before, that 

A SS SSf 

Now we introduce the relations Sf = Sfo + Sfm and S = S^ -\- Sms to treat 
separately the pieces that existed in the massless case from the innovations of 
the massive case. Thus 



\5Sq \SSms.~ SSfo , SSfm .„„^ 



Given the assumption that the new terms Sfm and Sms correspond, this equa- 
tion separates into the familiar part 

6Sfo _ A SSo . ^ 

and the new part 

^Sfm _ A 5Sms I ~ 
8^^ ~ 2 8fi^ 
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Using the gauge invariance result for the massless part of the free field action, 
etc. as above, one derives the form of 5*0 to be 



So = Si[gfi„,u] + S2, (33) 

with 5*2 as in the massless case. 

Assuming the free field mass term to be quadratic in the gravitational po- 
tential, one can assume the form 

Its variational derivative is 

Changing to the bimetric variables gives 



'-{[C + D]fj'X-Cg^VABV^-'-Dg^r^fj'^ + E6:fj'X-E9^rBrAy, 

(34) 

the term with is a reminder that the detailed coefficients depend on the 
space-time dimension, which has already been assumed to be 4. The new part 
of the action for the massive case satisfies 

SSfm _ A SSms I - 

so that becomes, using the explicit quadratic form of Sfm, 



^{[C + D + iE]f,>X - Cg^VABr" - D-g^fi^n'X - E~g^ fi'i.f,^^) = ^^Iff- 



Using det{f]^) = -y/^ , one finds the useful result 



dy/—rj 1 



dfj^^ 1 — Aw 



One now needs to try to find an appropriate form for Sms- It appears 
impossible to find anything to contribute a term of the form Cg^rjABV^" , one 
can show that this apparent impossibility is genuine, because the term is not 
an exact differential, whereas the remaining terms, with suitable coefficients, 
are. Thus C = 0. For the remaining terms, a natural form to try is Sms = 
J d^x{Pgl^fi'g + Q)y/--rj for mispecified real numbers P and Q to be determined 
shortly. One notes that a pure -^/^ piece that gives a cosmological constant 
plays no role here in this derivative, and can be already included in Sq. 

Equating ^ '^/-T |g with ^^f- and equating coefficients determines several 

i/j. I n 

of the constants. Equating the coefficients of the ^—^ffA terms gives Q = 
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+ 4i?)(l - Aw) (recalling that = 327rG'). Equating the coefficients 

327rG • 



of the \/—riVB^ASiy terms gives P = • Equating the coefficients of the 



—Tjfj'^'if^g^ terms gives P = ^^J^^q^' ■ Using all three results together gives 

E ^ 



Aw-r 
P= ^ 



AnG 

AnG 



Q---^^^ (35) 



Thus one has 

{Aw + 3)D D 

= 

and 

Combining the algebraic piece of Sq with S^s gives 

When the gravitational potential vanishes, 5*0(3 ought to vanish as well. Im- 
posing this condition to zeroth order gives B = Recalling that the 
goal is to find a massive generalization of Einstein's theory (with no effective 
cosmological constant), we require the first-order term in 7^ to vanish as well. 
It follows that A = 2L»(1 - Aw) and B = j^. Thus the sign of the formal cos- 
mological constant term depends on the density weight of the potential chosen 
initially. One also expects the quadratic part of the algebraic piece of the action 
Saig to agree with the free field mass term S/m- After a binomial expansion and 
some algebra, one sees that this is the case. One has 

Cat, = D^r, [itn^l^fi'k - + HOT, 

where HOT involves cubic and higher order terms. 

Making a weak-field expansion of the full massive nonlinear action S to relate 

2 

the coefficient D to the mass m of the spin 2 gravitons shows that D — 
For non-tachyonic theories, one imposes D <0. A. helpful fact in this expansion 
is that using the field redefinition = V'^ — ^1^, the perturbation of the 
associated metric is, to first order, 

g^u = fiiJ.<^ - ^^liixv) + ■■■■ 

This choice of convention, inspired by ([11, 12]), makes the traceless part of 
the metric perturbation independent of the type of gravitational potential {e.g., 
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covariant vs. contravariant) as far as is readily possible. This goal explains 
some of the factors of 5 that appeared above. 

Combining all these results gives the total massive action S, which depends 

on the spin 2 graviton mass m and the density weight parameter w which 
controls the relative mass of the spin to that of the spin 2: 

S = j d'^Xy/^R{g) + SmatteriQuu, u] 

+^ J dH^Tgi-^ - H + V^l^w + 3] - V^g^^A), (37) 

for w ^ i. Note that the vanishing of the coefficient (1 — Aw) of the ^/—g term 
would be a catastrophe and so is forbidden; its forbiddenness appeared at the 
beginning of the derivation, where that value would yield a densitizcd metric 
describing only the null cone, not the whole metric tensor, and so would involve 
coupling only to the traceless part of the stress-energy tensor [64]. By contrast, 
the vanishing of the coefficient of the is permissible and implies conformal 
invariance. The first line in the action contributes the ingredients of General 
Relativity, the second line contributes nothing to the field equations (but some- 
thing to the Roscnfcld stress-energy tensor), and the third line contributes the 
mass term. These theories are all universally coupled, pace the claim [8, 9] that 
the Preund-Maheshwari-Schonberg theory has that feature uniquely. 

One readily sees that, on fixing the local Lorentz freedom, these theories 
fall within the 2- parameter Ogievetsky-Polubarinov family [11]. One needs to 
identify the parameter w with their — | and to notice that a cotetrad is roughly 
an inverse square root of the contravariant tensor g^'^ that they take as basic, so 
their parameter n here takes the value of — ^ . One can then use their calculation 
of the ratio of the spin mass mo to the spin 2 mass m: 



, -16w2-8w; + 3 
mo = m\ -—: : (38) 

(The calculation is mildly nontrivial because the density weight of the poten- 
tial affects the trace piece of the mass term.) Requiring the mass to be real 
(non-tachyonic) excludes a great many values, but leaves two intervals. (By 
contrast the densitized metric theories contribute a single interval [10].) The 
mass ratio is shown in the figure. One recalls that the ^—77 term may vanish, 
giving a massless spin and conformal invariance, but the v^— 5 term must not 
vanish. The permitted intervals are [— |, — 5] and [0, j). For w = — |, conformal 

invariance and masslessness for the spin obtain. For w = — ^'^^ ~ —.683013, 
mo = m. On this branch the spin mass rises with w; the value w = —■^ 
gives an infinitely massive spin 0, as indicated by the vertical asymptote in the 
graph, and reduces to the Pauli-Fierz pure spin 2 case to quadratic order. The 
case w = on the other branch is similar. The cases w = — i (weight h) and 
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Figure 1: Ratio of Spin Mass to Spin 2 Mass as Function of w; done with 
Mathematica. 

w = Q (weight 0) correspond, respectively, to the fourth and second pure spin 
2 hadron mass terms of Zumino [65, p. 492]. From the vertical asymptote the 
mass descends, the ratio reaching 1 dX w — ~ .183013. As w approaches 

i from the left, the spin mass goes to 0, but it never reaches because w ^ j. 
Higher values of w lead to imaginary masses. In contrast to the two metric- 
based famihes previously derived by universal coupling [10], the tetrad theories 
permit the spin mass to be heavier than the spin 2. Such flexibility expands 
the possible astrophysical and cosmological phenomenology [38,39] of massive 
gravities to theories derived using universal coupling. This flexibility might also 
be relevant to the stability issue, which is expected to depend on resonances 
[42] as appeared above; resonances in this case are determined by the ratio of 
the spin and spin 2 masses. 

Reconsidering the generalized Bianchi identities for both coordinate and lo- 
cal Lorentz freedom in terms of the bimetric variables, one can find a shortcut 
to the inference of a Lorenz-Lorentz-type auxiliary condition, that is, one sim- 
ilar to the electromagnetic condition d^A^^ = 0. Using both identities and the 
gravitational and matter field equations, one can infer that 

d,{fl^^\~9-w5^^f,t^\g)=i). (39) 

The massive theory's action splits into three parts, as noted above, which act 
differently here. The first part, 

j d'^Xy/^Rig) + Scatter [cit^u , u] , 

which contributes general relativistic terms to the field equations, contributes 
nothing here, because the fiat tetrad is absent. The second part 
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which contributes nothing to the field equations, contributes here a term with 
that is nonzero to {ri^j^\g — wS^fj^j^lg^ , but its divergence vanishes iden- 
tically. The third piece, the algebraic part, which gives the graviton masses, 

; J d^x{^{l - 4w) + ^[Aw + 3] - V^g^fj'X), 



8ttG . 

gives an interesting result: 

54.g,%-[l + H5p'^«) = 0. (40) 

This condition eliminates a spin 1 and a spin from the field content. It is closely 
analogous to the fact that the Lorenz-Lorentz condition dfj^A^^ = follows for 
massive Proca electromagnetism using the electromagnetic field equations and 
(to enforce charge conservation) the equations of motion for charged matter. 
Ogievetsky and Polubarinov [11,12], who did not employ a tetrad formalism 
(but rather a surprisingly capacious formalism that takes powers of the metric 
using binomial series expansions), derived massive gravities using such a spin 
limitation principle instead of universal coupling. 

The mass term is the only place where the tetrads, as opposed to the met- 
rics, appear essentially. Thus the mass term is the only piece for which the 
local Lorentz invariance identity gives interesting results. That identity can be 
written as 

6S ^ SS ^ 
9^lB—^+V^^[B—^=0. (41) 

Using the gravitational equations of motion gives 

^.[B^=0. (42) 

Inserting the contribution from the mass term gives the result 

~9p[Bfi% = 0, (43) 

or equivalently, 

9[^.Afi^] = 0; (44) 

the antisymmetrization applies only over the Greek indices in this equation. 
While the tetrads are asymmetric, expressions with either only Latin or only 
Greek indices are symmetric. 



5 Massless Universally Coupled Gravity with Den- 
sitized Tetrad 

I now turn to deriving first General Relativity, and then massive variants thereof, 
using the contravariant version of the derivation above. There is not a com- 
plete symmetry on account of the fact that a covariant Lagrangian density has 
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weight 1, not weight 0; thus one cannot simply swap indices up and down and 
interchange some + and — signs everywhere. But conceptuahy the process is 
identical, so the exposition will be more concise. 



5.1 Massless Free Field Action 

Let Sf be the action for a set of four free covcctor densities of density weight 
w, w J. The weight 2w densitized inverse metric tensor 7]^"^ is essentially the 
densitized tetrad squared: 

fj'"' = V^AV^'^v'h- (45) 

It is the inverse of r)^^ for the same value of w. 

Let us require that the free field action S / change only by a boundary term 
under the infinitesimal gauge transformation 

Sj'^ = -rAd.e + wrAd.C- (46) 
One obtains the identity 

There is also the local Lorentz invariance of the free theory: 

Using an arbitrary infinitesimal change of coordinates with compact support, 
one derives the generalized Bianchi identity for any invariant action. Using the 
matter and gravitational equations of motion, as well as the other (local Lorentz) 
Bianchi identity, one arrives at an expression for a conserved symmetric stress- 
energy tensor (density): 

* =°- 

One can impose universal coupling with the postulate 

S^_SSf_ A 6S 



One is free to make a change of variables from 7^ and f)^ to the 'bimetric' 
variables and fj^, where 



rA=rA+^7i. (51) 



Equating coefficients gives 



SS ,^ 6S ,^ SS 
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for 5?7^ and 

5S A SS 



for ^7^. The universal coupling postulate becomes 

A 8S SSf 

Multiplying by the flat tetrad, taking the divergence, and recalling the free 
field theory's Bianchi identity from gauge invariance gives 

SS s-u ~a I- 



d.{filj^\~g-wS;rAj^\~9]=0- (55) 



One can show that the resulting equation is, in terms of the non- weighted tetrad 



d^(vl^\9^=0. (56) 

If one assumes that curved tetrad g'^ appears only through the curved metric 
Qn^, as one knows occurs in Einstein's equations, then the flat tetrad also appears 
only through its corresponding metric. By the same reasoning as above, one 
obtains 

+ y^dW,^, (57) 

which gives Einstein's equations. 



6 Massive Tetrad Density Theories 

The tetrad derivation can be generalized in the presence of a mass term. Thus 
Sf — Sfo + Sfm and 5 = + Sms- Requiring 5/0 to change only by a boundary 
term under the variation 

Sl'X = -fi''Ad.^'' + wrAd.C- (58) 
One obtains the identity 
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There is also the local Lorentz invariance of the free theory: 



One can impose universal coupling with the postulate 

S^a' S^A^ 2SrA. ^ ' 

One is free to make a change of variables from 7^ and 77^ to the 'bimetric' 
variables g'^ and ry^. The universal coupling postulate becomes 

\ SS ,^ SSf 



Now one makes the corresponding splits Sj = Sjo + Sfm and S = So + S„ 
The massless parts Sfo and 5*0 go as in the massless case and yield terms that 
give Einstein's equations and terms that affect only the Rosenfeld total stress- 
energy tensor. It remains to consider the equation 



5Sfm _ A dS„ 



\9- 



Assuming the free field mass term to be quadratic in the gravitational potential, 
and recalling the exclusion of the C term above, one can assume the form 

Sfm = V^rAlBim^f,^ + Efl^f,^). 

Its variational derivative is 

Changing to the bimetric variables gives 

^ = ^-^{-Dnt + ^s'^Bt.nt - ^Ef^i + Eg^sV^vt)- (63) 
Thus one infers that 

^-^i-Df}^ + WB^lt - ^Ent + ErsV^f,^) = -^^15- 
One can show that 



~A. 



For Sms a natural form to try is Sms — J d'^x{Pggfi^ + Q)^J—r]. Employing 
this expression and equating coefficients of like terms determines several of the 
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constants. Equating the cocfScicnts of the \/—rifj^ terms gives Q = j^{D + 
iE){4:W — 1). Equating the coefficients of the y/—rj''lf^''lv9B terms gives 

AttG 

Equating the coefficients of the yZ—VV^V^ 9b terms gives 

E{l-4w) 
A-kG ■ 

It follows that E = Thus one has 

and 

Combining the algebraic piece of Sq with Sms gives 
A 



Salg = -TT^ / <fx^+2B / (fxyf^ 



8nG 



When the gravitational potential vanishes, 5*^;^ ought to vanish as well. Im- 
posing this condition to zeroth order gives B = ■ Recalling that the 
goal is to find a massive generalization of Einstein's theory (with no effective 
cosmological constant), we require the first-order term in 7^ to vanish as well. 
It follows that A = 2D{Aw — 1) and B = 0. One also expects the quadratic part 
of the algebraic piece of the action Salg to agree with the free field mass term 
Sfm- After a binomial expansion and some algebra, one sees that this is the 
case. One has 

Calg = (rAV^rBV^ + ^^^^ + HOT, 



Aw 

where HOT involves cubic and higher order terms. Making a weak-field expan- 
sion of the full massive nonlinear action S to relate the coefficient D to the mass 

2 

m of the spin 2 gravitons shows that D = For nontachyonic theories, one 

imposes D < 0. 

Combining all these results gives the total massive action S, which depends 
on the spin 2 graviton mass m and the density weight parameter w which 
controls the relative mass of the spin to that of the spin 2: 



^ WttG . 



SwG 



J d^V^iAw - 1) + v^[5 - Aw] - V^g'^fi^), (65) 
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Figure 2: Ratio of Spin Mass to Spin 2 Mass as Function of w; done with 
Mathematica. 

ior w ^ J. 

One readily sees that, on fixing the local Lorentz freedom, these theories also 
fall within the 2-parameter Ogievetsky-Polubarinov family [11]. One needs to 
identify the parameter w with their — | and to notice that a tetrad is roughly a 
square root of the contravariant tensor g'^'^ , so their parameter n here takes the 
value of ^. One can then use their calculation of the ratio of the spin mass 
niQ to the spin 2 mass m: 



/-16iy2 + 24w-5 

Requiring the mass to be real leaves two intervals. (By contrast the densitized 
inverse metric theories contribute a single interval [10].) The mass ratio is shown 
in the figure. The permitted intervals are (|, 5] and [1, |]. As w approaches j 
from the right, the spin mass goes to 0, but it never reaches because w ^ j. 

The spin mass then rises with w; for w — ~ .316987, mo — m. w — ^ 

gives an infinitely massive spin 0, as indicated by the vertical asymptote in 
the graph, and reduces to the Pauli-Fierz pure spin 2 case to quadratic order. 
The case w — 1 on the other branch is similar. The cases w = and w = 1 
correspond, respectively, to the third and first of Zumino's pure spin 2 hadron 
mass terms [65, pp. 490, 492]. From the vertical asymptote at it; = 1 the mass 
descends, the ratio reaching 1 at w = ~ 1.18301. For w = | conformal 

invariance and masslessness for the spin obtain. 

Reconsidering the generalized Bianchi identities for both coordinate and lo- 
cal Lorentz freedom in terms of the bimetric variables, one can readily find the 
Lorenz-Lorentz-type auxiliary condition. Using both identities and the gravita- 
tional and matter field equations, one can infer that 

9'^{v:.j^j9-n.S;r.j^Jg)=0. (67) 
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Only the mass term 

2 

^(v^(4«; - 1) + V^[5 - 4w] - V^-giXv^), 

gives an interesting result: 

9.m^ + [l-t«]5H^<5p = 0. (68) 

Because the mass term is the only place where the tetrads, as opposed to 
the metrics, appear essentially, it gives an interesting result using the equations 
of motion, namely, 

ff^'I^r?^! = 0, (69) 

or equivalently, 

-g^Xr^^ = 0. (70) 
Thus the potential is symmetric when expressed using indices of the same type. 



7 Conclusion 

Taking into account the results of the present paper and its predecessor [10], re- 
expressing old results using I = 2w as needed, one has these four one-parameter 
families (with isolated and non-isolated forbidden cases) of universally coupled 
massive gravities. For the covariant symmetric tensor (density) potential with 
weight —2w, one has the universally coupled mass terms [10] 

rn'^ 1 

.{^[l-Aw]+V^[Aw + l]--^g^.rn (71) 



with —j<w<j. For a contravariant symmetric tensor (density) potential 
with weight 2w, one has [10] 

-XV^i'iw - 1] - v^[4«; - 3] - - v^^^-fy^,), (72) 



ISttG^^ ' " " ' 2 

for i < w < |. As derived above, for the cotetrad (density) case with weight 
—w one finds 



SttG 



{^[l-Aw]+^[4w + 3]-^g^fj'X). (73) 



The permitted intervals for w are [— f,— 5] and [0, |). Finally, for the tetrad 
(density) case with weight w one obtains 



8^ 



(^[4«;-l] + y^[5-4«;]-v^5^j7;^). (74) 



The permitted intervals for w are ij,^] and [1, |]. All of them fit within the 
Ogievetsky-Polubarinov framework, after discarding the local Lorentz freedom 
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if necessary. The massive tetrad theories, unhke the massive metric theories, 
permit the spin mass mo to be heavier than the spin 2 mass m, not merely 
hghter or the same, and so have more flexible phenomenology. The greater 
the variety of universally coupled mass terms, the more probable it seems that 
some of them can escape dangerous resonances sufficiently to have a chance at 
stability. It is therefore important both to pursue maximal generality and to 
investigate the question of stability, whether analytically, numerically or both. 
In view of the long-range diflSculties experienced by General Relativity, such 
as dark matter and dark energy, there is also empirical motivation for careful 
exploration of the possibilities for massive gravity. 
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